Ward identities for disordered metals and superconductors 



Revaz Ramazashvili 

Materials Science Division, Argonne National Laboratory, Argonne, IL 60439, USA. 

Abstract 

This article revisits Ward identities for disordered interacting normal metals 
and superconductors. It offers a simple derivation based on gauge invariance 
and recasts the identities in a new form that allows easy analysis of the quasi- 
particle charge conservation (as e.g. in a normal metal) or non-conservation 
(as e.g. in a d-wave superconductor). 

I. INTRODUCTION 

Interplay of interaction and disorder remains one of the central topics in condensed 
matter physics. Given the complexity of the problem, constraints imposed by symmetries 
acquire particular importance. An example of such a constraint is given by Ward identi- 
ties. In the early days of many-body theory, Ward identities were used to establish key 
properties of the Fermi liquid theory. 1 In the context of the CPA approximation for a dis- 
ordered non-interacting metal, similar identities were derived by Velicky. 2 In the theory of 
superconductivity, Ward identities were used early on to establish gauge invariance of the 
electromagnetic response. 3 Subsequently, they were employed by D. Vollhardt and P. Wolfle 4 
in a self-consistent theory of the Anderson transition, by F. Wegner, 5 A. J. McKane and 
M. Stone 6 in the sigma model approach to localization, and by C. Castellani et al. 7 in an 
early treatment of an interacting disordered metal. Very recently, T. R. Kirkpatrick and D. 
Belitz 8 invoked the Ward identities in an attempt to resolve the issue of decoherence at zero 
temperature. 
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This paper revisits the Ward identities for disordered interacting normal metals and 
superconductors. Using gauge invariance, it derives the identities in a new form that makes 
quasiparticle charge conservation (as e.g. in a normal metal) or absence thereof (as e.g. in a 
d-wave superconductor) explicit. In a normal metal, the identity takes a particularly simple 
form: 



where Ara is the disorder average of the retarded-advanced charge density vertex correction 
at zero momentum transfer Q = p — p = and small frequency transfer Q = uj — uj' <C 
uj,uj', and S^,(w,p) is the imaginary part of the retarded quasiparticle self energy which is 
proportional to the quasiparticle scattering rate. The vertex Ara is closely related to the 
correlation function of the quasiparticle charge density, and the 1/ (uj — uj') behavior of the 
vertex at low frequency transfer and zero momentum transfer points to quasiparticle charge 
conservation and its diffusive propagation. By contrast, in a d-wave superconductor, the 
Ward identity reflects the fact that impurity scattering causes exchange of charge between 
the quasiparticle subsystem and the condensate, which leads to non-conservation of the 
quasiparticle charge. 

The structure of the paper is as follows. Section II gives a detailed derivation of the Ward 
identities for a disordered interacting normal metal. Section III briefly discusses the Ward 
identities for an s-wave superconductor in the approximation of a spatially uniform gap. 
Section IV derives the Ward identifies for a disordered d-wave superconductor in the same 
approximation, and Section V illustrates the meaning of the identity by explaining how, in 
a d-wave superconductor, the impurity scattering leads to exchange of charge between the 
quasiparticle subsystem and the condensate. Section VI presents a summary and a brief 
discussion of the results. 





cu — uj' 



II. WARD IDENTITIES FOR A NORMAL METAL 



Consider a disordered interacting normal metal with the Matsubara action 
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where ■0 + (■0) are the electron creation (annihilation) operators, A is the electromagnetic 
vector potential, 0(r, r) is the scalar potential and u(r) is the impurity potential. This action 
respects the continuous gauge symmetry 
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of which the sought Ward identities are a consequence. To establish the scheme used through- 
out the rest of this article, I present below a detailed derivation. 

Everywhere hereafter, only infinitesimal time dependent spatially uniform transforma- 
tions ip a (r,r) — > e lx ^ip a (r,T) will be considered. Under such a transformation, the Green 
function changes according to 

G afS (r, r', r-r')^ e^G Q/3 (r, r', r - r^e'^ 

and thus, to first order in x, its variation equals 

5G aP (r, r>, r - r') » i[ X (r) - x(r')]G a/3 (r, r', r - r'). 

On the other hand, the same transformation induces extra terms in the action due to the 
presence of the temporal derivative. Hence, the very same variation of the Green function 
can also be calculated by perturbation theory. The crucial point is that the four-fermion 
interaction term in the action is invariant under the gauge transformation and, therefore, 
does not contribute to the perturbative correction to the Green function. Thus, to first order 
in infinitesimal x, the same correction to G is equal to 

5G a/3 (x,x',T - t') = -i J dtdr(tl) a (x,T)^(r,t)^(r,t)^(x',r'))dtx(t). 

Equating the two expressions leads to the identity 

\x(r)-x(r')}G a p(x,x',T-r') = - J dtdr{^ a (x,T)^+(r,t)^(r,t)^{x! ,r'))d tX (t) 



for a given disorder configuration. Disorder averaging replaces the exact Green function on 
the left hand side by its translationally invariant average. The average on the right hand 
side can be presented as the product of the two average Green functions plus the vertex 
correction term and, for x — Xo etQT with f2 — > 0, the Fourier transformed identity takes the 
form 

G{iu + ifl,p) — G(iu,p) = iVLG{iu + iQ,p) [1 + A(iu,iu + iQ;p,p)] G(iu,p), 

where A(iu,iu + ifl;p,p) is the disorder average of the scalar vertex correction. At this 
point, two different types of identities can be derived: one for the retarded-advanced vertex 
correction A R a(u>,u + Q;p,p), and another one for the retarded-retarded vertex correction 
A RR (uJ,uJ + tt;p,p). 

A. The identity for the retarded-advanced (RA) vertex 

To obtain the identities for the retarded-advanced vertex, choose iu to be in the lower 
half-plane and iu+ifl in the upper half-plane. Then, upon analytic continuation iuj — > u±i0, 
G{iui) transforms into Ga(^ — i0), whereas G(iu + ifl) transforms into G R (u + Q + iO). The 
identity then takes the form 

G r \lu + n + i0,p)- G A \cu - i0,p) = n[l + A RA (co + Q,uj;p,p)]. 

The disorder averaged Green function reads G A ^ R (cu,p) = u — Y,a/r{u,p) — where 
T,a/r(uj,p) is the advanced/retarded self energy. Using the relation T,r(u,p) = Y, A (uj,p) and 
assuming that the derivative dJTiR/A{wiP) is non-singular, for small VL = u — uj' — > one 
finds 

-2iY!' R {u,p) = [u -w']Ara(w,uj';p,p). (1) 

Identifying 2T, R (u) with the scattering rate 1/r, one immediately recognizes in A R a the zero 
momentum transfer (Q = 0) limit of the charge density vertex D{uj — u', Q) 9 
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D(u - u', Q) = . 

i[uj — uj')t + DQ 2 r 

where D is the diffusion coefficient. For a non-interacting disordered metal, D(u — u/, Q) is 
commonly obtained by a direct calculation, 9 first finding self-consistently the impurity self 
energy, and then summing the ladder series for the vertex. In the presence of interactions, 
diagrammatic treatment becomes much more involved, while the present derivation appeals 
only to gauge invariance and is insensitive to turning on the interaction. 

B. The identity for the retarded-retarded (RR) vertex 

By contrast with the identity just derived, the identity for the retarded-retarded vertex 
can be found in textbooks, and I present its derivation here only for completeness. In this 
case, it is convenient to choose both iu and iuj + iQ in the same (say, the upper) half-plane. 
Upon analytic continuation and multiplication by G~ 1 {iu) and + iQ), the identity 

takes the form 

G R \u + tt + iO,p) -G R \u + iO,p) = + A RR (uj + n,uj;p,p)], 

which, to first order in Q — ■> 0, leads to the standard relation between the energy derivative 
of the retarded self energy and the retarded-retarded vertex 10 : 

<9 w £ r (cj) = A RR (u,uj;p,p). (2) 



III. WARD IDENTITY FOR AN S-WAVE SUPERCONDUCTOR 



In the Nambu notations, the BCS Hamiltonian of an s-wave superconductor reads 



H 



-J 



e -> 



T ^(P At 3 ) + 7i A(r) + r 3 e0 + t 3 u 



Here ^ = is the Nambu spinor, Tj are the Nambu matrices, and r denotes the 

center of mass coordinate of a Cooper pair. The pair field A(r) has been chosen real for the 
sake of simplicity. 



A gauge transformation takes the form \1> — >■ exp ir^x ^ an d, in addition to the stan- 

— # 

dard change of potentials A and <f>, has to be accompanied by the pair field transformation 
A — > A exp 2ij^X ■ One then proceeds the same way as for a normal metal, with two im- 
portant points to note. The first point amounts to the approximation of a spatially uniform 
gap which, along with the frequency, gets renormalized by disorder. The second point stems 
from the Nambu matrix structure of the theory: the vertices, that appear after disorder 
averaging of the perturbative expression for the Green function variation, are defined in 
the Nambu space and thus carry a Nambu index. For instance, the disorder averaged term 
arising from the temporal derivative of x — Xo ex P[^ r ] has the form 

J dy^{x)i>{y)iQ.T^{y)i>{x')) - f dydzdz'G(x, z)Xl [r 3 5(z - y)5(y - z') + (r 3 ) (z, y, z')} G(z', 

where the vertex correction (t 3 )(z, y, z') appears as a result of disorder dressing of the corre- 
sponding bare vertex, the latter being simply the Nambu matrix r 3 . The disorder averaged 
Green function has the form 1 

G{p,u) = [iujT - £(jo)r 3 - Atx] -1 , 

where u and A are the renormalized frequency and the gap amplitude, which yields the 
Ward identity 

[iuj - iuj'] T 3 — i A UJ + A w / r 2 = [iu - iu'] [r 3 + (r 3 )] - 2?A[r 2 + (r 2 )]. 

This, in turn, indicates a diffusion pole in the quasiparticle charge density vertex correction 
(73) jm : 

-2iY!' R (u,p) 

r 3 = {T3/RA, 



uo — UO' 



where Y!' R (uj,p) is the imaginary part of the retarded self energy renormalization of the 
frequency, the notation is chosen to coincide with the normal metal limit. 
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IV. WARD IDENTITY FOR A D-WAVE SUPERCONDUCTOR: 



In a d-wave superconductor, the situation turns out to be quite different. The BCS 
Hamiltonian of a d-wave superconductor reads 



where the pair field A(R, r) has been chosen real and having d-wave angular dependence on 
the relative coordinate r, and R denotes the center of mass coordinate of a Cooper pair. 

As in the s-wave case, the Hamiltonian respects the gauge symmetry, and the identities 
can be obtained similarly, with one important difference: because of the d-wave symmetry of 
the gap and its oscillating angular dependence, the gap amplitude A p , although suppressed 
by impurities, does not acquire a frequency dependent renormalization. Hence the disorder 
average of the quasiparticle Green function is 



Another important point is that the angular dependence of the gap leads to the appearance 
of a vertex correction (A P T2) on the right hand side of the Ward identity, which assumes the 
form 



As for an s-wave superconductor, S^(o;,p) is the retarded self-energy renormalization of the 
frequency: uj = uj — E. Due to the d-wave symmetry of the gap and its oscillatory angular 
dependence, (A p t 2 )ra oc (t 3 )ra, which leads one to conclude that the vertex correction 
( t 3)ra has to remain finite as uj — uj' — > 0. Hence, in a disordered d-wave superconductor, 
the quasiparticle charge is not conserved. Note that, upon transition to the normal state, 
the quasiparticle charge diffusion mode re-appears, as can be seen seen by sending A v to 
zero in Eq. (3) and identifying (t 3 )ra with Ara((v,u';p, p) of Section II. 



H = 




G(iuj,p) = [iuj -TiAp- r 3 £ p ] 




(3) 
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V. QUALITATIVE ARGUMENT FOR A D-WAVE SUPERCONDUCTOR 



The absence of quasiparticle charge conservation in a d-wave superconductor can be 
understood based on a simple argument going back to the studies of charge imbalance 
relaxation in superconductors. 11 I reproduce the argument here for the sake of completeness. 
Consider the Bogolyubov quasiparticle creation operator: 



Impurity scattering is elastic, i.e. it conserves the quasiparticle energy E p = + . In 

an s-wave superconductor with uniform gap, A p is a constant and, in the absence of the 
Andreev scattering that turns £ p into — £ p , the energy conservation implies conservation of u p 
and v p . Hence the impurity scattering conserves the particle-hole content of a quasiparticle, 
and this leads to the effective charge conservation - even though a Bogolyubov quasiparticle, 
being a superposition of a particle and a hole, does not have a well defined charge quantum 
number. The same conclusion can be reached by considering directly the expectation value 
of the quasiparticle charge Q p : 



In an isotropic s-wave superconductor, the gap does not vary around the Fermi surface and 
hence, in the absence of the Andreev processes, Q p is conserved by the impurity scattering, 
which leads to the charge diffusion pole. 

By contrast, in a d-wave superconductor, the gap A p is strongly anisotropic. Thus, 
even in the absence of the Andreev scattering processes, neither Q p nor the moduli of 
the Bogolyubov factors u p and v p are conserved: impurity scattering changes the particle- 
hole content of a quasiparticle. Physically, this means that the impurity scattering induces 
exchange of charge between the quasiparticle subsystem and the condensate. 

Indeed, this quasiparticle charge non-conservation is not a consequence of the d-wave 
symmetry of the gap, but rather of the gap anisotropy around the Fermi surface, and is 
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present not only in other superconductors with non-trivial symmetry, but even in s-wave 
superconductors with anisotropic gap. However, in the latter case, the effect is small in the 
measure of the relative gap anisotropy, which is itself reduced by disorder. As a result, the 
quasiparticle charge non-conservation appears only at time scales that are long compared 
with the scattering time. By contrast, in a d-wave superconductor, the gap anisotropy is 
large, and the quasiparticle charge changes at the time scale of order the impurity scattering 
time, which eliminates quasiparticle charge conservation at any time scale beyond the elastic 
scattering time. 

VI. SUMMARY AND DISCUSSION 

In this article, I revisited the Ward identities for superconductors and disordered inter- 
acting normal metals, and presented a simple derivation based solely on gauge invariance. 
The identities were recast in a new form that made quasiparticle charge conservation (as in 
a normal metal or an isotropic s-wave superconductor) or absence thereof (as in a d-wave 
superconductor) explicit. Using the Ward identities, I showed how, in a d-wave supercon- 
ductor, impurity scattering causes exchange of charge between the quasiparticle subsystem 
and the condensate, thus leading to the quasiparticle charge non-conservation. 

Transparency of the Ward identities is particularly appealing in comparison with mi- 
croscopic approaches. The simplicity of the identities is insensitive to the strength of the 
impurity potential or to whether disorder has to be treated in the Born or in the unitary 
limit - or to the presence of interaction. By contrast, to achieve a controllable approxima- 
tion even in the Born limit, microscopic calculations, e.g. for a d-wave superconductor, have 
to resort to rather complex methods and/or unrealistic approximations, such as expansion 
in the inverse number of gap nodes. 12 

It is my pleasure to thank A. J. Leggett, M. R. Norman and especially Yu. M. Galperin 
for helpful discussions. Work supported by the U. S. Department of Energy, Division of 
Basic Energy Science-Material Science under contract No. W-31-109-ENG-38. 
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